Current Rectification and Seebeck Coefficient of Serially Coupled Double Quantum 
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The transport properties of serially coupled quantum dots (SCQDs) embedded in a matrix con- 
nected to metallic electrodes are theoretically studied in the linear and nonlinear regimes. The 
current rectification and negative differential conductance of SCQDs under the Pauli spin block- 
ade condition are attributed to the combination of bias-direction dependent probability weight and 
off-resonant energy levels yielded by the applied bias across the junctions. We observe the spin- 
polarization current rectification under the Zeeman effect. The maximum spin-polarization current 
occurs in the forward bias regime. Such behavior is different from the charge current rectification. 
Finally, the Seebeck coefficient (5')of SCQDs is calculated and analyzed in the cases without and 
with electron phonon interactions. The application of SCQDs as a temperature detector is dis- 
cussed on the basis of the nonlinear behavior of S with respect to temperature difference across the 
junction. 
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1. Introduction 

Serially coupled quantum dots (SCQDs) exhibit the 
transport properties of current rectification due to the 
Pauli spin blockade, negative differential conductance 
(NDC), nonthermal broadening of the tunneling cur- 
rent, and coherent tunneling in the Coulomb block- 
ade regime. 1-3 ^ Although many theoretical works have 
been devoted to investigating these phenomena, they 
still can not explain the transport properties of SCQDs 
systematically. 4 " 8 ' Sun e al calculated the tunneling cur- 
rent of SCQDs in the Pauli spin blockade using the 
Keldysh-Green function technique. 4 ' The procedure in- 
troduced in ref [4] to solve high-order Green functions 
arising from electron Coulomb interactions can not re- 
solve the quantum paths of SCQDs. In Rcfs. 5-8, the 
master equation was used to calculate the tunneling cur- 
rent of SCQDs. However, in these works, cases are re- 
stricted to t c -C T, where t c and T denote, respectively, 
the interdot hopping strength and tunneling rate between 
the electrodes and the quantum dots (QDs). 

Here, a closed-form expression for the tunneling cur- 
rent of SCQDs with a finite interdot hopping strength 
enables the analysis of current rectification arising from 
coherent tunneling with a spin blockade and the NDC of 
tunneling current resulting from off-resonant energy lev- 
els. The effect of Zeeman energy splitting on tunneling 
current is investigated to clarify the behavior of spin- 
polarization current. In addition, the Seebeck coefficient 
(S) of SCQDs is calculated in the cases without and with 
electron phonon interactions (EPIs). In the absence of 
EPIs, we propose how to use SCQDs as a temperature 
detector on the basis of the nonlinear Seebeck coefficient. 
When the SCQDs are embedded in a phonon cavity 9-11 ', 
it is possible to manipulate EPIs to control the electrical 
conductance and Seebeck coefficient of junction systems. 

2. Formalism 

Because we consider nanoscale semiconductor QDs, 
the energy level separation of between QDs is much larger 
than their on-site Coulomb interactions and thermal en- 
ergies. One energy level for each quantum dot is con- 



sidered in this study. The two-level Anderson model in- 
cluding EPIs is employed to simulate the SCQD junction 
system shown in the inset of Fig. 1(a). The Hamiltonian 
of an SCQD junction 12 ' is given by H = Hq+Hcqd+Ht 



H = $^e*4,cr *,* + ^2 e k b ia b k.« (1) 

k.(T k.(T 

where the first two terms describe the free electron gas of 
the left and right metallic electrodes. a\ a (b\ J) creates 
an electron with momentum k and spin a with energy 
€k in the left (right) metallic electrode. Vkj, (i = 1,2) 
describes the coupling between the metallic electrodes 
and the first (second) QD. d\ a {dt t „) creates (destroys) 
an electron in the £-th dot. 



Hcqd = ^2 E e n LlT + ^2 Uem,, 



(2) 



where Et is the spin-independent QD energy level and 
n>£,<7 — &\ a di.o- Notations Ug and Utj describe the in- 
tradot and interdot Coulomb interactions, respectively. 
ttj describes the electron interdot hopping. describes 
the EPIs 



H T = ojqJc + ^2 Vtfrii^^ + c), 
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where ujq is the phonon cavity frequency and is the 
coupling strength of EPIs. A canonical transforma- 
tion can be carried out to remove on-site EPIs, that is, 
H new = e s 'He s , where S = -J2 La Q e n e ^ - c). 13 < 14 ) 
In the new Hamiltonian, we have the following effec- 



tive physical parameters: Vi 
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V k ^ c *-<=>, Ef = E t - X 2 e cj , Uf = U t - 2A|w , 
Ufa = U t ,j - 2AiA 2 w , t e ej = t^e-^- x ^- c \ and 
\e = Cli/u>o. Under the canonical transformation, the 
coupling strengths between the electrodes and the dots, 
on-site energy levels, intradot Coulomb interactions, in- 
terdot Coulomb interactions, and electron interdot hop- 
ping strengths are renormalized by EPIs. If we consider 
a special case of Ai = — A2 = A, we have 

H T = w c f c + to{ni,<r - n 2 ,a){c 1 + c). (4) 

a 

This special case of Eq. (4) was already considered in refs 
15 and 16. For the case of Eq. (4), we have an effective 
electron interdot Coulomb interaction Ui.j+2X 2 ujo, which 
is always repulsive and enhanced with increasing EPIs. 16 ) 

To decouple the EPIs of H newi we take the mean-field 
average to remove the phonon field arising from — c, 
which is < ea;p A * (ct ~ c) >= exv -\*i™tb?\u 1 {2k B T p )\ _ 0n 
the basis of such a mean-field average, we see a reduc- 
tion of v fei£ exp-5 A ? cot ' l2 K/(2fe B T p )] = v k jX t and inter- 
dot hopping strength t c exp-^ ( - Xl -^ 2coth2 ^^ (2kBT ^ = 
t c Xej. This leads to the redefinition of the renormalized 
resonance width of each energy level of the QD. T p is the 
phonon temperature. 

Using the Keldysh-Green function technique and ne- 
glecting the phonon -assisted tunneling process, 12 > the 
tunneling current of an SCQD is given by 

J =f / deT(e)[f L (e)-f R (e)], (5) 

where T(e) = (7i2(e) + 72i(e))/2 is the transmission fac- 
tor. /L = i (j R =2 )(e) = l/^-^Ms))/^^*) + 1] denotes 
the Fermi distribution function for the left (right) elec- 
trode. The left (right) chemical potential is given by 
Ml (Mi?)- Ml — i-iR — eAV a , where AV a denotes the ap- 
plied bias. Notation T L r R \ denotes the equilibrium tem- 
perature of the left (right) electrode, e and h denote 
the electron charge and Planck's constant, respectively. 
Tej(e) denotes the transmission coefficient, which can be 
calculated by the on-site retarded Green function and 
the lesser Green function. The transmission coefficient 
has the following expression: 

* r|(e)r^' m (e) 

where Im means taking the imaginary part of the function 
that follows, and 

GJ imi<r (c)=p m /( W -S m ). (7) 

r l=L(i),ii(2)( e ) = r t x h where r ^ = E fc Vg tt 5(e - c fc ) de- 
notes the tunnel rate from the left electrode to dot A 
(E\) and from the right electrode to dot B (E 2 ), which 
is assumed to be energy- and bias-independent for sim- 
plicity, in = e — Ef + iTg/2. We can assign the follow- 
ing physical meaning to Eq. (6). The sum in Eq. (6) 



is over eight possible configurations labeled by to. We 
consider an electron (of spin a) entering level i, which 
can be either occupied (with probability Ni yS ) or empty 
(with probability 1 — N^). For each case, the electron 
can hop to level j, which can be empty (with probability 
dj = 1—Nj^—Nj^+Cj), singly occupied in a spin a state 
(with probability b^ s — N^ s — Cj) or spin a state (with 
probability b^ a — Nj >tr — Cj), or in a double-occupied 
state (with probability Cj). Thus, the probability fac- 
tors associated with the eight configurations appearing 
in Eq. (6) arc p 1 = (1 - N e ^)aj, p 2 = (1 - N tt9 )b jt9 , 
P3 = (1 - Ni^)bj >a , Pi = (1 - N t ^)cj, p 5 = Ne^cij, 
Pe = N£^b 3 ^, p 7 = Ni^bj >a , and p s = N^Cj. E m in the 
denominator of Eq. (7) denotes the self-energy correc- 
tion due to Coulomb interactions and coupling with level 
j (which couples with the other electrode) in configura- 
tion to. We have E x = t 2 /nj, E 2 = + 1 2 /(^ - Uf), 
E 3 = U^+e/i^-Ulz), E 4 = 2Ul j +i 2 /(» j -U?-Ul e ), 
S5 = Ui + t 2 /( P] - E 6 = Uf + Uf. + t 2 /( P] - 

U* - U? e ), E 7 - Ut + Ut d + t 2 /(^ - 2U? e ), and 
E 8 = U c t + 2Ufj + t 2 /(^ - Uf - 2U? e ). Note that 
t = t% = t c X ej . Here r^ m = -2ImEj denotes the 
effective tunneling rate from level I to the other elec- 
trode through level j in configuration m. For example, 
Tf = ~2lmt 2 /^ = t 2 Ty[{e - E]) 2 + (Tp2) 2 }. It is 
noted that T^' m has the numerator Tj for all configu- 
rations. Furthermore, G r ea (e) = X) m =i Ge, m ,cr( e ) is sim " 
ply the on-site single-particle retarded Green function for 
level I as given in Eq. (A16) in Ref. 12, and G r tma (e) 
corresponds to its partial Green function in configuration 

TO. 

The probability factors of Eq. (7) are determined by 
the thermally averaged one-particle occupation number 
and two-particle correlation functions, which can be ob- 
tained by solving the on-site lesser Green functions: 

fde ^ TUe(e) + T e j ' m f j (e) T 

N i,v = - / — 2^ pe + pe.m ^G^^e), 

(8) 

and 

Ce = ~ / V 2^ r e + r e,m ImG^^e). (9) 

Note that I ^ j in Eqs. (6), (8) and (9), which are valid 
under the condition of t c /Ue <C 1. We will study the 
transport properties of SCQDs on the basis of Eqs. (5), 
(8), and (9). 

3. Results and discussion 

3.1 Current rectification 

To numerically calculate the tunneling current of SC- 
QDs without EPIs (Ai = A2 = 0), we adopt the intradot 
Coulomb interactions of Ug = U and tunneling rates of 
T L = T R = T. We consider these conditions of ho- 
mogenous intradot electron Coulomb interactions, and 
symmetrical tunneling rates for simplicity. All energy 
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scales are considered in units of Tq. On the other hand, 
771(2)6 A V^, is employed to describe the energy shift aris- 
ing from the applied bias AV a across the junction. That 
means that Eg is replaced by eg = Eg + r]geAV a , assum- 
ing the right electrode is grounded. On the basis of the 
experiment in Rcf. 1, we adopt rji = 0.6 and rj2 = 0.4. 
Although the factor rjg depends on the QD shape, ma- 
terial dielectric constant, and location, we assume that 
rji is determined by the QD location, that is r/g = Lg/L, 
where Lg is the distance between the grounded electrode 
and the £th QD, and L is the separation distance between 
the left electrode and the right electrode. 

We plot the tunneling current of SCQDs under the 
Pauli spin blockade condition (Ex + £/ 12 = E 2 + U2) 
shown in the inset of Fig. 1(a) for three different in- 
terdot Coulomb interactions at fc^T = 1 r , t c = 0.1 r , 
U = 30 r , and T = T . The three curves corre- 
spond to (a) U12 = and E\ = Ef, (b) U12 = 5 T 
and E 1 = E F - 5 T , and (c) U 12 = 10 T and 
Ei = Ef — 10 r . Ef denotes the Fermi energy of the 
electrodes. The current rectification and negative differ- 
ential conductance (NDC) of SCQDs are observed. The 
maximum tunneling current is suppressed in the pres- 
ence of interdot Coulomb interactions. We note that the 
Jmax,Rl 'Jmax,F ratios are near 2. J ma x,R and J ma x,F 
are the maximum current in the backward bias and for- 
ward bias, respectively. J max ,Rl 'J ma x,F = 2 is in very 
good agreement with an experimental observation. 1 The 
results shown in Fig. 1 indicate that the Pauli spin block- 
ade condition is not attributed to interdot Coulomb in- 
teractions, but to intradot Coulomb interactions. On 
the basis of Eq. (6), the Pauli spin blockade resonant 
channel of Ei + U12 = E 2 + U2 is determined from the 
probability weights of p 2 = (1 — Ni^)(N 2id — c 2 ) and 
P5 = N 2 ^(l — Ni, s — Ni i<r + Ci), which result from T12 
and 72i, respectively. Figure 1(b) shows the occupa- 
tion number Nu 2 ^^ = -^1(2), a = -^1(2) an d two parti- 
cle correlation functions Cu 2 ) for Ui 2 = 10T . The P2 
of 7i2 and p 5 of 721 are also plotted as indicated by 
the dotted line and dashed line in Fig. 1(b), respec- 
tively. Because J max ,F and J ma x,R occur, respectively, 
at eAV a = 3r and eAV a = -4.5r , we have p 2 = 0.28 
and P5 = 0.3 at eAV a = 3Tq and P2 = 0.24 and p§ = 0.65 
at eAV a = — 4.5r . Their sum is P2+P5 = 0.58 at J ma x,F 
and P2 + Ph = 0.89 at J ma x,R- This demonstrates why 

J max 

, R is larger than J ma x,F- In the reversed bias, the 
increase in p§ with respect to eAV a results from the en- 
hancement of (1 — Ni i9 — Ni i(7 + ci) (dot A empty), which 
provides an increased probability of tunneling of electrons 
in the right electrode. This result of p$ approaching one 
in the high backward bias indicates that the current spec- 
tra shown in Fig. 1(a) are determined not only by the 
probability weights of P2 and P5 but also by off-resonant 
energy levels. We find that the off-resonant energy level 
is a key reason to observe the NDC behavior of SCQDs. 
As a consequence, the QD energy level broadening will 
significantly influence the maximum currents in the for- 
ward and reversed biases. 3 ) It is worth noting that the 



phonon- assisted tunneling process arising from EPIs can 
not be ignored in the high bias regime. 7 ) Some current 
structures of SCQDs 1 ) in the high-bias regime can be well 
explained by the phonon-assisted tunneling process. 7 -* 

For the applications of SCQDs in spintronics, it is cru- 
cial to measure the spin configuration of each electron 
in individual QDs. SCQDs have a functionality of spin- 
charge conversion. 1 ~ 3 ' However, it is not easy to measure 
the small magnitude of tunneling current in the weak 
interdot coupling limit of t c /T <C 1. Although tunnel- 
ing current can be enhanced by increasing t c , we con- 
sider how the behavior of current rectification is influ- 
enced for t c /T fts 1. To clarify the above question, we 
plot the tunneling current of SCQDs at various t c val- 
ues in Figs. 2(a)-2(c). Other physical parameters are the 
same as those for the curve with U12 = 10r shown 
in Fig. 1(a). The maximum currents labeled J ma x,F 
and Jmax.R are shifted toward a higher bias when t c 
increases. The J ma x,R/ Jmax,F ratio slightly decreases. 
The three maximum currents are J moXi ir=10.4, 167.8, 
and 328, which correspond to eAV a =3, 3.6, and 4.8 To 
, respectively.(J maXj fj=21.6, 319, and 599 correspond to 
-4.5, -5.6, and -7.4Fo, respectively). The dashed lines 
shown in Figs. 2(a)-2(c) show the contributions aris- 
ing from only the resonant channels II2 of T12 and II5 
of 72i- These dashed lines very close the solid lines at 
the small applied bias (between J ma x,F and J m ax,R) but 
not at the large applied bias. From the results shown 
in Figs. 2(a)-2(c), the contributions of II2 and II5 still 
dominate the trend of current spectra. These two res- 
onant channels n 2 and II5 have the two poles e± = 
E F +0-5eAV a +iT/2±0.5y/{0.2eAV a ) 2 + At 2 c . Their bias- 
dependent probability weights p 2 = (1 — Ni^ s )(N 2 _ s — c 2 ) 
and P5 = A^2,ct(1 — Ni t9 — Ni, a + ci) are determined by 
the occupation numbers shown in Figs. 2(d)-2(f). The 
effect of interdot hopping on Ni and N2 is enhanced with 
increasing t c . The results shown in Fig. 2 indicate that 
the current rectification behavior of the SCQD is not de- 
stroyed under the condition of t c /T w 1. 

Although the SCQD system has the functionality of 
spin filters under the Pauli spin blockade condition, the 
tunneling currents in Figs. 1 and 2 do not exhibit 
the spin-polarization current. The spin-polarization cur- 
rent of SCQDs was theoretically studied in Ref. 4 in 
which the authors considered the spin-bias and ferro- 
magnetic electrodes. Our study introduces the Zeeman 
effect arising from a local magnetic field to yield the 
spin-polarization current under the Pauli spin blockade 
condition of Ei + Ui 2 = E 2 + U 2 . Here the energy 
level of each QD depends on electron spin. That is 
,Ei,a- — Eg ±giHBBa = Eg =p Az t e t<T is considered, where 
ge denotes the g-factor, \ib is the Bohr magneton, and 
B is the magnetic field. For simplicity, we assume the 
homogenous g factor gi = g < 0. We plot the spin- 
dependent charge current J a (J -a) and spin-polarization 
current J p = J a — J_ CT in Fig. 3. Other physical param- 
eters are the same as those for the curve of Ui 2 = 10r 
shown in Fig. 1. The increase in spin-polarization cur- 
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rent is observed with increasing Zeeman splitting. Unlike 
charge current, the maximum spin-polarization current 
in the forward bias is larger than that in the backward 
bias. That is J p , m ax,F > J P ,max,B- Figure 3(c) shows 
the spin-dependent occupation number which determines 
the spin-dependent probability weights. So far, we have 
considered Tp — T R = T in Figs. 1-3, where the charge 
currents are generated only by the applied bias. To study 
the thermoelectric effect of SCQDs, we will consider the 
case of T L ^ T R to investigate the Secbeck coefficient. 
3.2 Seebeck coefficient 

If the SCQD junction system is in an open circuit, 
an electrochemical potential will form in response to a 
temperature difference across junction [see Eq. (5)]; this 
electrochemical potential is known as the Seebeck volt- 
age (Secbeck effect). The Seebeck coefficient (amount 
voltage generated per unit temperature gradient ) is de- 
fined as S = A/// AT, where Ap = pp — PR = eAV a 
and AT = T L — T R are the voltage difference and tem- 
perature difference across the junction, respectively. Re- 
cently, many studies have been devoted to the investiga- 
tion of the thermoelectric effects of SCQDs in the linear 
response regime. 17 " 18 ) Because the efficiency of thermo- 
electric devices can be enhanced with a large temperature 
difference, 19 ^ it is important to clarify the behavior of 
the Seebeck coefficient at a large temperature difference 
AT. Previous studies focused on the nonlinear thermo- 
electric properties of individual QD systems, 20 > 21 ) which 
may not be readily realized from the experimental point 
of view. The thermal resistivity of the SCQD junction 
system can be larger than that of a single QD system. 
This feature allows the SCQD system to maintain a rel- 
atively large temperature difference across the junction. 
In this section, we study the Seebeck coefficient in the lin- 
ear and nonlinear regimes. An analytical solution of S in 
the linear response regime gives useful guidelines for un- 
derstanding the behavior of S in the nonlinear response 
regime. 

In the linear response regime, Eq. (5) can be rewritten 
as J = L eAV a + TifceAT. The thermoelectric coeffi- 
cient L n is given by 

2g r E \ n 

Ln = ~h J den€ hcoshi((e-E F F )/(2k B T)y (10) 

where the transmission coefficient T(c) is calculated 
under the equilibrium condition. We define S(T) = 
AV a /AT = (fc B /e)(-Li/ J Lo). Figure 4 shows the Sec- 
beck coefficient (S) of SCQDs with identical QDs as 
a function of temperature at various QD energy levels 
Eq — Ep = A = 20,30,40,50 r . The negative S indi- 
cates that the diffusion electrons pass through the reso- 
nant channels oi E ±t c , E + Ui2±t c , E n + U + [7 12 ±t c , 
and Eq + U + 2U\2- These resonant channels result from 
the four configurations of pi, p^, pe, and p$ in Eq. (7). 
The maximum S(T ) increases with increasing A. T 
denotes the temperature at which S has a maximum 
value. In addition, T is shifted toward to the high- 
temperature regime. The dashed lines calculated using 



Sq(T) = — A/T are employed to fit these solid lines. 
These fitting curves show good agreement with the solid 
lines in the high -temperature regime. The characteristic 
of So = —A/T, which is independent of the tunneling 
rates, was also determined in the case of a single QD 
junction. 20 -* In the appendix, we derive the formula of 
So(T) = —A/T for the noninteraction case when the 
conditions of ^/fc^T <C 1 and r/fc^T <C 1 are satisfied. 
Such results imply that the behavior of So(T) is not sensi- 
tive to electron Coulomb interactions. This is attributed 
to the very small contribution of resonant channels in- 
volving electron Coulomb interactions that are far above 
Ep. Note that the behavior of S(T) becomes complicated 
if when Eg — Ef = A < 0, because S{T) will involve 
many physical parameters such as electron Coulomb in- 
teractions, the electron interdot hopping strength t c , and 
tunneling rate. 17 - 1 

To calculate the Seebeck coefficient in the nonlinear 
response regime, we numerically solve Eq. (5) by consid- 
ering the condition of J = 0. Figure 5 shows the elec- 
trochemical potential and Seebeck coefficient as a func- 
tion of temperature difference AT for different detuning 
energies A at fc^T^ = 40 To. The curve with trian- 
gles(considering A = 50 Tq) neglects the shift of the 
QD energy level yielded by the electrochemical poten- 
tial. That is, r/i = r)2 = 0. Ap, increases with increasing 
AT. When AT > 0, Ap = eAV a is negative. On the 
other hand, Ap is positive when AT < 0. By Compar- 
ing the blue solid line with the curve with triangles in 
Fig. 5(a), we find that under the forward temperature 
bias (Tl > Tr), the case of rp = 7/2 = needs a larger 
electrochemical potential (Ap) to balance diffusion car- 
rier flow mainly through the resonant channels of Eq ± t c 
from the left electrode to the right electrode. This is be- 
cause the resonant channels are always kept in the case 
of rji = 7] 2 = 0. Once r/i ^ 0, the electrochemical po- 
tential will generate the off-resonant channels and the 
carrier diffusion flow will be blocked. Under the reverse 
temperature bias (UbTr > ksTp = 40r ), the difference 
between the blue solid line and the curve with triangles 
is small even though AT = — 20Fo. This indicates that 
diffusion carrier flow does increase too much with increas- 
ing temperature bias. Figure 5(b) shows the Seebeck co- 
efficient corresponding to the electrochemical potentials 
shown in Fig. 5(a). The behavior shown in Fig. 5(b) 
can be roughly described by 5* = S /{1 - 0.5AT/T). If 
AT/2T <C 1, we have T R = T L -Ap/S . This result may 
be useful for the application of temperature detectors. 22 ) 

According to the definition of S = (fcs/e)(— Li/L ), S 
is related to the electrical conductance G e = eL . We 
further examined the relationship between these ther- 
moelectric response functions. Figure 6 shows the elec- 
trical conductance (G e ) of SCQDs with identical QDs 
(Eg = Eq = Ep + 30Tq — eV g ) as a function of gate volt- 
age for different values of EPIs (fie) at a low temperature 
ksT = ir . Note that we have T p = T in the linear re- 
sponse. The coupling strength of EPIs can be tuned by 
the phonon cavity, which modulates the phonon density 
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of states to change A^ = f^/wo. 9 n ) In the absence of 
Xi = 0, there are eight resonant channels of E ± t c , 
E a + U 12 ±t c , E a + U + U 12 ±t Cl and E a + U + 2U 12 ± t c , 
which are labeled from eVgi to eV g g. The first peak (eV g \) 
is shifted toward a low gate voltage with increasing A, 
which is attributed to the shift of the QD energy level 
E\ r = En — \ 2 (jJ - Because of the reduction of the in- 
terdot Coulomb interaction Uf - — Uej — 2Xi\ 2 ujo, there 
are only six peaks and four peaks in the Ai = X 2 = 0.5 
and Ai = A 2 = 0.7. For Ai = A 2 = 0.7, Uf^ almost 
vanishes. Therefore, four peaks correspond to Ef ± t c , 
and El + Uf ± t c . For the homogenous coupling of EPIs 
Ai = A2, t% = t c is independent of A. The separation 
between the peaks corresponding to the bonding and an- 
tibonding states is invariant. 

In Fig. 6, we considered the case of Ai > and \ 2 > 0. 
It is possible to have a positive Ai and a negative A 2 by 
considering QDs at particular locations in the phonon 
cavity. For instance, dot A is out of phase from dot B 
under phonon fields [see Eq. (4)]. This will lead to a re- 
duction of intradot Coulomb interaction and an enhance- 
ment of interdot electron Coulomb interaction. Figure 7 
shows the electrical conductance (G e ) and Seebeck coef- 
ficient (So) as a function of gate voltage for different A 
values. For Ai = 0.5 and \ 2 = —0.5, the spectra of G e 
and So are significantly changed owing to the reduction of 
interdot hopping strength [i2 e -(A 1 -A 2 ) 2 coth 2 (o;o/(2fe B T))] 
and the enhancement of interdot Coulomb interactions 
Ui.j — 2AiA 2 o;o. The Seebeck coefficients shown in Figs. 
7(c) and 7(d) correspond to Figs. 7(a) and 7(b), re- 
spectively. The Seebeck coefficient shows a change in 
sign arising from the bipolar effect which is electron-hole 
asymmetrical, where holes are defined as the empty states 
below the Fermi energy of electrodes. 17 -* 

4. Summary and Conclusions 

Intradot Coulomb interactions are more important 
than interdot Coulomb interactions in the observation 
of current rectification spectra in the Pauli spin block- 
ade. The current rectification is not restricted in the 
weak interdot hopping limit (t c /T <C I). From the ex- 
perimental point of view, the release from the weak inter- 
dot hopping restriction is very useful for measuring the 
tunneling current of SCQDs. Unlike the charge current 
rectification, the maximum spin-polarization current is 
obtained in the forward bias regime. The Seebeck coef- 
ficients in the linear and nonlinear response cases were 
studied. The universal behavior of the Seebeck coeffi- 
cient (So = — k B A/eT) may be useful for the application 
of temperature detectors. EPIs provide an extra degree 
of freedom to control carrier transportation. 
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Appendix A: Seebeck coefficient 

The tunneling current given by Eq. (5) can be analyt- 
ically calculated by contour integration. For simplicity, 
this appendix only considers the case without electron 
Coulomb interactions. 
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t 2 c [fL(e)-f R (e)} 



\( e - ei +iT L /2)(e-e 2 + iT R /2)-t^- 

(AI) 

fL(R)( e ) is the Fermi distribution function of the left 
(right) electrode, ti = Eg + rjieAV a . In the linear re- 
sponse regime, we take eAV a — > and /c^AT — > in Eq. 
(AI) and we have J = L eAV a + Lik B AT, where 



L n /a = / 

J —1 



dy 



y n /cos 2 (y/2) 



(y - - uX)(y - uj-)(y - cj* ) ' 

(A2) 

w± = (E±-E F )/(k B T), and E± = 



Where a = ih(k B T)* 

e + iV/2 ± Eg. We define e = (E 1 +E 2 )/2 = E , T = 
(T L + T R )/2 = To, and Eg = ^((E x - E 2 )/2) 2 + if. The 
Seebeck coefficient is defined as So(T) = —Li/L (kB/e) 
for considering the condition of J = 0. 

By using contour integration, L and L\ can be eval- 
uated in terms of polygamma functions. 23 ) Because we 
are interested in the case of fcsT/r 3> f and A = 
(Eq — Ep) » t c , we have the thermoelectric response 
functions L and L\ 



2e ii 



7rr /2 



h k B Ttl + (To/2) 2 cosh 2 (A/2k B T) 



(A3) 



and 



Li 



2e i 2 



7rr /2 



h k B T 2 1 2 + (r /2) 2 cosh 2 (A/2k B T) ' 



(A4) 



On the basis of Eqs. (A3) and (A4), we obtain the 
electrical conductance (G e ) and Seebeck coefficient (So) 
with the following expressions 



G e - 



and 



2e 2 t 2 



irTo/2 



h k B T t 2 c + (Tq/2) 2 cosh 2 (A/2k B T) 



(A5) 
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*o = (A6) 

The Seebeck coefficient of Eq. (A6) depends on the de- 
tuning energy and equilibrium temperature. 
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Figures and Figure Captions 




eAV a /r 



FIG. 1: (a) Tunneling current as a function of applied bias 
for various interdot Coulomb interactions (U12) at tempera- 
ture UbT = ITq. The tunneling currents are in the units of 
Jo = 2eTo/h. (b) Occupation number (Ne.) and two -particle 
correlation function (ce) in the case of U12 = WTq. 
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FIG. 2: Tunneling current as a function of applied bias. 
Diagrams (a)-(c) show the interdot hopping strengths t c = 
0.1,0.5, and 1 To, respectively. Diagrams (d)-(f) showing N 
and ce correspond to diagrams (a)-(c), respectively. Note that 
ci is not plotted in diagrams (d)-(f) because of the very small 

Cl. 




FIG. 3: Spin-dependent tunneling current as a function of 
applied bias: (a) A z = 0.3 T and (b) A z = 0.7 r . A z = 
gfiBB/2. The spin-polarization currents (dashed lines) are 
defined as J p — J a — J- a . (c) Spin-dependent occupation 
number as a function of applied bias at Az = 0.7Fo. 




FIG. 4: Seebeck coefficient (So) as a function of temperature 
for different detuning energies (A = 20,30,40,50 To) in the 
linear response regime. The physical parameters Ui — U — 
30ro, Uij = 5To, r_L = Tr = To, and t c = 3Fo are adopted. 
Dashed lines are calculated using Sf = —A/T, which is given 
by Eq. (A6). 




FIG. 5: (a) Electrochemical potential Afj, and (b) Seebeck 
coefficient (S) as a function of temperature difference AT for 
different detuning energies at temperature of the left electrode 
Tl ~ 40Fo. Other physical parameters are the same as those 
in Fig. 4. 
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FIG. 6: Electrical conductance (G e ) as a function of applied 
gate voltage for different electron-phonon interactions (EPIs) 
at k B T = ir , E e = 30r -eV g , U e = 60r , and U t ,j = 10r . 
Diagrams (a)-(c) show the results for Ai = A2=0, 0.5, and 
0.7, respectively 




FIG. 7: (a)-(b) Electrical conductance (G e ) and (c)-(d) See- 
beck coefficient (So) as a function of applied gate voltage for 
different EPIs at k B T = ir . Solid lines: Ai = A 2 = 0, 
dashed lines: Ai = — A2 = 0.5. Other physical parameters are 
the same as those in Fig. 6. 



